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Abstract. Measurement of heavy quark production is one of the tools used to 
investigate the matter produced in extremely hot and dense conditions in heavy 
ion collisions at RHIC. The PHENIX experiment has measured mid-rapidity 
transverse momentum spectra of electrons. After subtracting the photonic 
background contribution, the electron spectra are mainly due to semileptonic 
' decays of hadrons containing heavy quarks and therefore provide a measure- 

\^ , ment of heavy quark production and its energy loss in hot and dense matter. 

■ This paper will present the technique used by the PHENIX experiment and 

recent results on heavy quark production in p+p , d+Au and Au+Au collisions 
at ^/sWn =200 GeV. 



1. Introduction 



o 

^ ■ Keywords: heavy quark, charm, bottom, RHIC, PHENIX 

Ph! PACS: 24.85.+P, 25.75.-9, 25.75.Dw 

> 

X, 

^ , Heavy quarks can be produced by quark-antiquark interaction forming a virtual 

gluon decaying into a cc or 66 pair. But the most important production mechanism 
is gluon fusion. Therefore, their spectrum is sensitive to initial parton distribu- 
tion. The heavy quark can bind to a lighter quark and form open charm (D) or 
bottom (B). Due to the heavy mass of the reaction products, a large momentum is 
transfered during the process, corresponding to a short distance interaction leading 
to a small coupling constant in the QCD framework. Their cross section can 
be calculated in QCD in terms of ag^" (pQCD). Fixed Order plus Next-to- Leading 
Logarithms calculations (FONLL) [ 1 , with c —t D and b B non-perturbative 
fragmentation functions, have shown reasonable agreement with Tevatron bottom [ 
121 13] and charm [j4j data at ^/s — 1.8 TeV and midrapidity. Since heavy quarks are 
likely produced during early stages of a nucleus- nucleus collision and interact with 
the media, they can be used as a powerful probe for nuclear and deconfinement 
effects. 
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Nuclear modifications are generally accounted for by a constant a in cr = ajsrA". 
There was no observation of any non- linear A dependence (a ^ 1) in open charm 
measurements at fixed target E769 (250 GeV tt+A) [0, WA82 (340 GcV tt+A) 
[El and E789 (800 GeV p+A) [Q experiments. However, WA78 (320 GeV tt+A 
in xp > 0.2, or small x region) has shown a significant non- linear target de- 
pendence {a — 0.81 ± 0.05) in D total cross section. Observations of charmonium 
suppression in p+A collisions [|^ ^1 and PHENIX d+Au data [^J are relevant 
to the question of shadowing of the small x gluon distribution in the nucleus [I13|. 

Some enhancement at high px of hadron production relative to p+p collisions 
(a > 1) has been observed in SPS [ El Ej and RHIC experiments [ E] some 
enhancement at hight . This pT broadening is known as Gronin effect [ I17j and it 
is attributed to the multiple scattering of the incident partons. 

A strong high pT suppression has been observed in PHENIX Au-|-Au collisions 
at y^SNN =200 GeV at midrapidity for 7r° and charged particles [ 1181 1191 1^ , 
despite the small enhancement in d-f-Au collisions measured in the same apparatus 
[ 1161 I21j . This behavior is well described by energy loss via gluon radiation in a 
hot and dense matter ES]. The gluon radiation occurs for angles 9 around the 
parton incidence higher than 0c = m/ E. Therefore, heavy quarks should have less 
energy loss due to the dead cone around the incident parton [I24|. 

Particle anisotropy is largely used to study the thermal conditions of the initial 
stage in non-central nucleus- nucleus collisions. In particular, the second component 
v2 of its Fourier expansion, called elliptic flow, is sensitive to early pressure. Posi- 
tive v2 is a indication of the particle thermalization [ 25 . A heavy quark anisotropy 
measurement can test the mass dependence of the elliptic flow which is an indica- 
tion of the media density. 

The amount of charm remaining after the hadronization stage in A+A collisions 
is an important information to constrain the weight of open charm recombination 
contribution in J/ip spectra. Direct measurement of D mesons at RHIG in d+Au 
collisions with px range up to 3 GeV/c was reported by STAR GoUaboration [I26j. 

In this paper results from electron decays of heavy quarks obtained in p+p , 
d-|-Au and Au-|-Au collisions at ^/snn— 200 GeV measured at mid-rapidity range 
in PHENIX Experiment will be shown. Most of the systematics is dedicated to 
extracting the photonic contribution from electron spectra by using two methods. 
Gomparisons of the methods are presented for consistency considerations. 

The data sets used in this paper can be seen at Table 
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2. Electron Measurement and Extraction of Heavy Quark Con- 
tribution. 

2.1. PHENIX Apparatus and Electron Identification. 

PHENIX Detector 




West Beam View East 

Figure 1: PHENIX central arm apparatus used for electron identification. 

The PHENIX detector has two central arms spectrometers (Fignj covering a 
pseudo-rapidity \ri\ < 0.35 and azimuthal angle A(p = 2 x tt/2. The total geometri- 
cal acceptance is 35%. A very important advantage in PHENIX central arms is the 
small amount of material. The total radiation length is Xq — 1.1% with Multiplic- 
ity Vertex Detector (MVD), during Runl-3, and Xq = 0.4% without MVD, from 
Run4 until now. 

The collision is triggered by a Beam Beam Counter (BBC, 3.1 < It^I < 3.9) and 



Table 1: Integrated luminosity used in each electron analysis at PHENIX. 



Run 


Projectile 


Integrated 






Luminosity 


Run3 


P+P 


350 nb-^ 


Run3 


d+Au 


2.74 nb-^ 


Run4 


Au-I-Au 


241 ^lb-^ 
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Zero Degree Calorimeter (ZDC, I77I > 4) coincidence in 92+3 q% of Au+Au inelastic 
cross section. The multiplicity measured by the BBC and energy from neutrons 
detected in the ZDC provide the centrality measurement of the collision. Tracks 
for charged particles with momentum higher than 0.2 GeV/c are recognized by the 
Drift Chamber (DC) and the pad chamber (PCI) with resolution Ap/p ^ 1%. Only 
tracks matching the Electromagnetic Calorimeter (EmCal) clusters within 2 apos 
are used. 

For the electron identification, we require a Cherenkov ring pattern around the 
track projection on the Ring Cherenkov Detector (RICH) and an energy deposited 
in EmCal cluster E corresponding to the particle momentum ^2^^ ' since hadrons 
deposit only a fraction of their energy in the calorimeter. The remaining contam- 
ination from particles with momentum lower than the RICH pion threshold (4.9 
GeV/c) comes from hadrons sharing Cherenkov rings with an electron. This contri- 
bution is statistically subtracted by using miscorrelated tracks (association of tracks 
with swapped Z coordinate in RICH). 

The electron spectra are corrected taking into account geometrical acceptance 
and reconstruction efficiency. They are calculated in a detailed GEANT-based [ 
I27| simulation previously tuned with fully reconstructed beam pipe conversion elec- 
trons [ED and the same dead areas as during data acquisition. The total electron 
acceptance times efficiency is almost uniform at 14% for momentum above 1 GeV/c. 
The multiplicity effect in detection efficiency is estimated by embedding simulated 
electron raw signal into real data. The fraction of reconstructed simulated electrons 
corresponds to approximately 20% for the most central events and 2% for periph- 
eral. These numbers are taken into account in the final efficiency. 

The final inclusive electron spectra for p+p and Au+Au data is shown in figure 



2.2. Cocktail Method for Photonic Contribution Estimation. 
The relevant contributions for the electron spectra are : 
• photonic source 

— Dalitz decays of light neutral mesons or internal conversion : 

— conversion of photons from collision point in material (beam pipe, MVD, 
air) : 

* light neutral mesons, i.e. : tt^ ^ 77 ^ 7e+e~ 

* direct photons 

— virtual photons: 7* e+e~ 
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• weak kaon decay K^z (/v * —t ir^e^Ve) 

• dielectron decays of vector mesons : p, oj, — s- e^e" 

• thermal radiation 

• semileptonic decays of D, and B 



The most important contribution comes from tt". The function E ^3^° is pa- 

rameterizated to tt" and charged pion spectra measured at PHENIX (Fig. [2J|. 

E is an input for a GEANT based hadron decay generator including the 

PHE&IX material for conversion estimation. 

TT° source function 



T^lnif 9.63/30 

C 403.2 ±92.91 

pO 0.7039 ± 0.05088 

a 0.3185 ± 0.05518 

b 0.07479 ± 0.04245 

n 8.318 ± 0.2076 
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6EANT based Monte Carlo for PHENIX 
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Figure 2: Parameterization of tt'^ spectrum function from tt'^ [ 18^ and charged 
hadron measurements [I^Oj at PHENIX. 



Assuming tot scaling and h/ir'^ ratio constant at high px, we use E 



d' a 



1^ 



P 



for 



other light hadrons h by replacing pT 
spectrum by the ratios : 



Pt + to| — TO^o and normalizing the 



• Tj/n^ = 0.45 ±0.10 
, ,//7rO = 0.25 ±0.08 

• Lu/n" ^ 1.0 ±0.3 

• p/tt° = 1.0 ±0.3 
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• (/./ttO = 0.40 ±0.12 

The spectrum obtained in this method agree whitin 20% for all pr range in rj 
and (f> yield measured at PHENIX . Direct photons [EH] and kaon measurements 
[1301 at PHENIX are used in the 7 and kaon decay generator. All electron source 
cocktail result are shown in figure |31 




(a) p+p data f I31| (b) Au+Au preliminary data 



Figure 3: Inclusive electron spectrum and relevant electron sources. 



2. 3. Converter Method. 

The conversion method is based on the installation of additional material for a short 
period of each Run for photonic source estimation. The "electron converter" is a 
thin brass tube (1.7% Xg) surrounding the beam pipe at 29cm. 

The electron yields without (ivf°""-°"*) and with (ArC"""-*") the converter 
from photonic (iVJ) and non-photonic (iV"°"^''') source are : 
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Figure 4: Data from Run2 [ ISE] (a) yields with (open circles) and without 
(filled circles) converter, (b) Ratio of converter in/out (Rcn) yields and ratio of 
photonic yields (R-y) determined by simulation, (c) Ratio of photonic/non-photonic 
yields (Rnp) and kaon contribution. 



During data acquisition with converter, e « 2.1% of electrons are lost. The 
converter contribution factor R^y is estimated for each photonic source by a full 
GEANT based simulation with and without the converter. The simulation agrees 
whithin 2% when compared to fully reconstructed conversion electrons in real data. 
The parameterizatcd function E "^Jlp" (Fig. ^ is used for the input in this sim- 
ulation. Likewise the cocktail method, uit scaling and tj/tt'^ ratio are applied to 
obtain the rj input. For the other mesons we assigned R!^ — since their total 
contribution correspond to only 6% of the photonic source at 3 GeV/c and their 
Dalitz branching ratio are similar to rj. The total R^ is a sum of the photonic 
sources weighted according to the hadron ratios. 

Figure 0] shows the momentum dependence of converter/non-converter yield 
(Rcn) and photonic/non-photonic yield (Rnp)- The non-photonic yield contains 
contribution from weak kaon decay and vector mesons. These contributions are 
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subtracted in the same manner as for the cocktail method. A good agreement is 
obtained between these different methods and data sets (Fig. [3)). 



Non-photonic component | 
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Figure 5: Photonic subtraction method and Run comparisons for Au+Au non- 
photonic electron spectra. 



3. Heavy Quark Results. 

In this section only electron spectra and yields from semileptonic decays of heavy 
quarks are reported. These data are usually mentioned as non-photonic electron, 
since most of the subtracted electron contribution are internal and external 7 con- 
versions. But in the data shown here, small contributions from kaons and vector 
mesons are also removed, according to the method described in the previous section. 



3.1. p+p data and pQCD results. 

The heavy quark production obtained in p+p collisions is reported in [ 31 . In 
figure El the spectrum is compared to leading order PYTHIA and FONLL pQCD 
calculations. 

The phenomenological parton shower program PYTHIA 6.205 used a modified 
set of parameters from central and bottom hadroproduction [ and CTEQ5L 
parton distribution functions [133- In order to accomodate NLO contributions, we 
use a scale factor K=3.5 in leading order contribution. The calculation describes 
well the data up to 1.5 GeV/c. 

Above 2 GeV/c the electron spectrum suggests contributions beyond FONLL 
pQCD calculation [El- Higher order terms and/or jet fragmentation can play a 



Electrons and Heavy Quark at PHENIX Detector 



9 



di 

10" 

E 

"a. 10'^ 
■a 

10"* 
10^ 

10"' 
10"* 

i ' 

s 

6 
5 
4 
3 
2 



°0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 

[GeV/c] 

Figure 6: (a) Non-photonic electron spectrum measured in p+p collisions with 
PYTHIA and FONLL [Elj predictions, (b) Data and FONLL ratio including error 
band and PYTHIA comparison. 

role in higher momentum range. 

The total cross section of charm quark-antiquark pair production was estimated 
by using the electron spectrum shape the PYTHIA extrapolation. The PYTHIA 
function including charm and bottom contribution was fitted to data with px > 0.6 
GeV/c. Extrapolating the charm contribution to pr = GeV/c, the mid-rapidity 
charm cross section obtained was ddce/dy = 0.20 ± 0.03(stat) ± O.ll(sys) mb. The 
rapidity integrated cross section correspond to CTcc = 0.92 ± 0.15(siai) ± 0.54(s?/s) 
mb. The same estimation, performed by FONLL, gives cTcc = 0.256lQ^4g mb. The 
systematics of these calculations is described in 



3.2. "Cold" Nuclear Effects in d+Au results. 

There is no indication of nuclear effect in d+Au collisions, since its preliminary 
non- photonic data scale with Ncoii (Fig. 0). The same behavior was observed in 
STAR fully identified D meson production [ I26| and other experiments covering 
large x region (see sectional. 
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(a) PRELIMINARY 60% < centrality < 88% (b) PRELIMINARY 40% < centrality < 60% 



(e +e-)/2 

Non-photonic 

Systematic error 

PHENIX ppfit 




^d+Au (0-20%) (e"+e-)/2 

-~- Non-photonic 

Systematic error 

PHENIX ppfit 




(c) PRELIMINARY 20% < centrality < 40% (d) PRELIMINARY 0% < centrality < 20% 

Figure 7: Preliminary non-photonic electron invariant cross section per binary col- 
lision in d+Au for different centralities and p+p best fit [^^. 



3.3. Energy Loss in hot and dense matter formed in Au+Au collisions. 

The heavy quark yield at mid-rapidity range scales with Ncoii (Fig. |Hla). This 
observation can be considered as an experimental verification of the binary scaling 
of a point- like pQCD [ Kffij . However, medium effects can influence the momentum 
distribution. 



Clear suppression is observed in central events at high px when Au-f Au data 
is compared to nuclear overlap integral < Taa > scaled p+p data (Fig. |Hlb). The 
quantification of this effect is obtained by Raa = (Tf^)xda"+ calculations (Fig. 
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N„„ Pt [GeV/c] 

(a) (b) 



Figure 8: Non-photonic electron : (a) Integrated yield (0.8 GeV/c < momentum 
< 4.0 GeV/c) versus Ncoii [EH [EH- (b) Run4 preliminary Au+Au spectra for 
different centralities compared to nuclear overlap integral < Taa > scaled p+p data. 



For px >2 GeV/c the heavy quark Raa is surprisingly similar to that for tt" 
observed by PHENIX [EIj since the dead cone effect is neglected for lighter quark 
made mesons. It is commonly considered the time-average transport coefficient q 
denoting the average squared transverse momentum transferred from a hard parton 
per unit path length while traversing the medium. In the figure 1^ the spectrum in 
most central events and the charm energy loss expected for three different q [ I87| 
are presented. 

Studies performed by S. Wicks et al. [ i38j including charm-t-bottom contribu- 
tion and assuming gluon density dNg/dy = 1000 are present in figure 1^ In one of 
the scenarios the elastic energy loss is considered in the final state spectra besides 
the gluon radiation found in light meson jet quenching. The suppression does not 
seem to be affected by the expected significant contribution from bottom at mo- 
mentum above 3 GeV/c. 



3.4. Elliptic Flow. 

In electron anisotropy analysis the photonic contribution is estimated by the con- 
verter method fsection f2.3|l for momentum below 1 GeV/c and the cocktail method 
for momentum above 1 GeV/c, based on v2 from 7r° PHENIX data. The prelimi- 
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Figure 9: Non-photonic electron Raa for different centralities. In panel (e) are 
shown theoretical expectations from charm (1-3) [ I37| and charm-|-bottom (radiation 
only and radiation-|-elastic bands) contributions [I38|. 



nary results and comparison with quark coalescence model [ I39| with/without charm 
flow are shown in figure [TUI This plot is consistent with published elliptic flow [I40|. 
but with more statistics. In the charm px range collective motion is observed as 
for light quark made mesons. The observed charm flow is consistent with partonic 
level thermalization and high density of the matter formed. The decreasing of v2 
at high pt can be explained by a reduced flow from bottom dominance in this range. 



4. Outlook and Expectations. 

PHENIX detector measured heavy quark production for p+p , d+Au and Au-fAu 
at mid-rapidity range for 0.8 < < 5 GeV/c by using its electron decay chan- 
nels. Photonic contribution to the electron spectra was calculated with two different 
methods which are in good agreement. 
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Figure 10: Non-photonic electron elliptic flow data and comparison with quark 
coalescence model with/ without charm flow. 



The next to leading order pQCD and PYTHIA calculations agree with p+p data 
at low pt and underestimates it at high px- The total charm- ant icharm cross section 
is (Tec = 0.92 ± 0.15{stat) ± 0.54(s?/s) mb and was estimated from a well parameter- 
izated PYTHIA extrapolation to p+p data. 

The total non-photonic production scales with the number of binary collisions. 
Medium effects were only observed in Au-|-Au collisions where the number of non- 
photonic electrons per number of collisions is suppressed at high pT- The intensity 
of this suppression is comparable to what was observed for lighter mesons. This 
result represents a puzzle for the mass dependent gluon radiation scenario which 
explains reasonably well the energy loss in tt" and charged particle jet quenching 
observed at RHIC. The bottom contribution does not seem to affect the suppression. 

Measurement of v2 points to a charm flow supporting a very high density in 
the early stages of the collisions. The reduced v2 at high momentum, despite the 
poor statistics in this range, can be explained by the contribution of a more massive 
particle, perhaps bottom. 

Run4 RHIC luminosity allows us to reach a more extended px region. Data 
from high luminosity Cu-|-Cu, p+p , and lower energy ^/snn = 62.4 GeV collisions 
are coming soon. They will constrain production and medium effects models. The 
future installation of a vertex detector will provide a better signal/background non- 
photonic electron measurement and it will make possible to disentangle charm from 
bottom production, since B has a significant larger lifetime than open charm. 
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